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SPHERICAL  PROPAGATION  MODELS  FOR  MULTIPATH-PROPAGATION  PREDICTIONS 

by 

L.P.  Ligthart 

Delft  University  of  Technology 
Dept,  of  Electrical  Engineering 
Mekelweg  4 
2628  CD  Delft 


Summary 

VA 

Multipath  fading  leads  co  a limitation  in  the  availability  and/or  reliability  of  microwave  links.  To 
study  the  propagation  mechanism  under  fading  conditions  propagation  models. based  on  ray  theory, above  a 
spherical  earth  have  been  developed  and  compared  to  the  well-known  planar  propagation  model  above  a "flat- 
tened earth".  The  reasons  for  studying  spherical  propagation  models  are  a)  to  avoid  the  limitation  of 
small  elevation  angles  in  planar  propagation  models,  b)  to  investigate  the  computed  and  measured  path 
delay  on  the  microwave  line  of  sight  links,  and  c)  to  set  up  an  analytical  approach  for  the  spherical  pro- 
pagation model  without  numerical  difficulties  and  cc  -puting-time  intensive  procedures. 

In  this  paper  computational  results  of  the  models  are  shown,  including  an  accuracy  analysis,  and  the 
use  of  spherical  propagation  models  is  illustrated  for  surface  duct  layers  above  water. 

rs 

1 • Introduction 


Multipath  propagation  can  cause  large  variations  in  the  received  signal  level.  i‘o  prevent  multipath 
fading  on  terrestrial  microwave  radio  links  two  diversity  techniques  are  frequently  used: 

- frequency-diversity, 

- height-diversity. 

Another  possibility  for  multipath  fading  reduction  is  derived  from  the  angle  of  arrival  dependency  of 
incoming  rays  at  the  receiver  antenna.  This  third  diversity  technique  looks  promising  when  use  is  made  of 
electronically  controlled  limited  scan  antennas  (Ref.  1).  Hybrid  reflector  antennas  for  multipath  fading 
reduction  have  determined  the  German  and  the  Dutch  contribution  to  the  COST  204  project  (European  Co-Opex- 
ation  in  the  field  of  Scientific  and  Technical  Research-Commission  of  the  European  Communities) . The  pro- 
ject is  entitled  "Phased  Array  Antennas  and  their  Novel  Applications”  (Ref.  2). 

In  The  Netherlands  the  Dr.  Neher  Laboratory  of  the  Dutch  PTT  and  the  Microwave  Laboratory  of  the  Delft 
University  of  Technology  have  collaborated  on  a research  program  to  investigate  the  possibilities  of  using 
limited  scan  antennas  on  digital  microwave  line  of  sight  links.  To  come  up  with  antenna  specifications  for 
fading  reduction,  group  delay  requirements  in  two  ray  propagation  models  have  been  considered.  Firstly, 
the  planar  propagation  model  has  been  used.  This  model  approximates  small  elevation  angles  and  can  be 
found  by  an  "earth  flattening"  coordinate  transformation  to  preserve  relative  curvature  between  rays  and 
the  earth.  Because  accurate  time  delay  information  is  needed  for  specifying  the  wanted  antenna  character- 
istics spherical  models  have  also  been  investigated. 

In  Section  2 some  attention  is  paid  to  ray  parameters  in  planar  and  spherical  models,  asstiming  a con- 
stant refractive  index  gradient  over  a given  height  segment. 

In  Section  3 results  cf  time  delay  computations  for  both  models  and  of  the  accuracy  analysis  for  the 
spherical  model  are  shown. 

In  Section  4 a modified  logarithmic  refractive  index  profile  is  assumed  in  order  to  study  the  pro- 
pagation mechanism  through  surface  ducts  above  water  under  different  spherical  model  conditions. 


2.  Ray  parameters  for  planar  and  spherical  models 


In  the  model  with  planar  earth  (Figure  1)  the  vertical  refractive  index  profile  m is  assumed  to  be 
dependent  on  z only  and  is  characterized  by  the  transformation 


dn(r)  ~ dm(z) 
dr  dz 


_1_ 

R 


1 


B" _ r 

a 


where 


n = refractive  index  in  the  spherical  model  with  earth  radiua  and  with  assumed  radial 


dependency  r 


Ml  H li  li 
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n = refractive  index  in  the  planar  model 

B"  = refractive  index  gradient  per  meter  in  the  planar  model 
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Figure  1.  Flanar  propagation  model  Figure  2.  Ray  parameters  in  planar  propagation  model 


Segmentation  is  carried  out  in  such  a way  that  B"  = constant  per  segment  p (Figure  2) . The  ray  para- 
meters are  functions  of  the  refractive  index  profile  and  6.  In  the  planar  model  the  ray  parameters  per 
segment  become 


ft 


E 


E 


i 
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= height  difference  of  the  ray  = segment  thickness 
= horizontal  distance  of  the  ray  along  the  earth  surface 
= path  length  of  the  ray 
= path  delay  of  the  ray 


If  a ray  reaches  a maximum  or  minimum  height  within  the  segment  an  extra  segment  is  introduced  based 
on  these  heights.  The  whole  radio  path  is  described  by  summing  the  ray  parameters.  The  appeal  of  the  planar 


In  the  model  with  spherical  earth  (Figure  3)  the  refractive  index  profile  is  characterized  by  the 
transformation 


dn(r)  J_+_L_b 

dr  dr  R R, 

f a f 


1 


where  N = refractive  index  in  the  spherical  model  with  arbitrary  fictive  radius  Rf 
B'  = refractive  index  gradient  per  meter  in  the  spherical  model 


Segmentation  is  illustrated  in  Figure  4 where  B1  = constant  per  segment  p.  The  ray  parameters  in  the 
sDherical  model  become 


- h = segment  thickness 

P sph 

- d = R,  . ti  = horizontal  distance  of  the  ray  along  the  earth  surface 

p h f P J 

- Sp  p = path  length  of  the  ray 


m-1  ■ ■ 


lift-  u rta 


Different  analytical  approximations  for  the  ray  parameters  have  been  derived  for  B'  > -N(rf)/rf, 
B’  < -N(r^)/rf  and  for  the  rays  which  reach  maximum  heights.  Per  segment  analytical  expressions  can  be 
obtained  by  making  approximations  in  the  ray  theory,  especially  of  the  term 


Using  Snell's  law  for  a spherical  layered  medium  the  ray  path  constant  C becomes 

C = N(rf)  . rf  . sin9 

For  B'  < -N(r^)/r^  the  term  I yields 

I « V^(r  2 - r 2)K  1 
m f p 

r is  the  maximum  r_  value  the  ray  would  reach  if  the  refractive  index  gradient  B'  = constant  were  extended 
m f p 

to  above.  This  means  r > r^.  The  factor  K is  chosen  so  that  for  r > r0 
m2  p m2 

2 2 2 2 2 
(r  - r9  ) K » r.V(r9)  - CT 
m 2 p 2 2 


and  for  r = r„ 
m 2 

K = -N2 (r. ) - B'  . N(r . ) . r, 
p 1 p 1 l 

For  B'  > -N'rf)/rf  the  term  I yields 

I s /k'  (r  2 - r2)  + K" 
p f 1 p 

K'  and  K"  are  chosen  so  that  the  approximation  is  optimal  for  r,  = r;  K'  and  K"  become 

p p * f 1 p p 

K'  = N2  (r. ) + B'  . N(r, ) . r. 

p 1 p 1 1 

2 2 2 
K"  = tr(r.)  . r,  - <T 

p 1 1 

3 . Path  delay- computations  and  accuracy  analysis 

To  make  comparisons  between  planar  and  spherical  models  computer  programs  have  been  developed.  The 
input  parameters  in  the  programs  are 

- link  geometry  (i.e.  antenna  heights,  distance  between  antennas,  earth  radius) 

- refractive  index  profiles  linearized  in  a given  number  of  segments 

- elevation  angles  of  transmitted  rays 

For  an  incoming  ray  at  the  aperture  of  the  receiver  antenna  the  results  are 

- path  length  and  path  delay  T 

- single  of  arrival  a 

- amplitude  and  phase 

In  case  the  elevation  angles  of  transmitted  and  received  rays  are  within  1 degree  the  main  differences 
between  the  models  are  mostly  found  in  path  delay  and  amplitude  computations.  Therefore,  in  this  .ection 
results  of  path  delay  computations  are  compared  for  the  spherical  model  with  Rf  = Ra  and  the  planar  model. 
For  simplicity  a constant  refractive  index  gradient  B'  and  thus  a constant  B"  has  been  assumed  over  the 
atmospheric  height  region  of  interest.  For  the  spherical  model  this  height  region  is  subdivided  into  :;eg- 


144 


ments  with  a constant  thickness  of  10  meters.  The  antenna  heights  at  transmitter  and  receiver  equals 
100  m above  the  earth.  Three  different  B*  (and  B")  values  are  considered: 


a.  B1  = -5  * 10-7  [m-1]  ■*  B"  = -3.4  * 10-7  [m-1] 

a a 

b.  B£  = -.78  * 10-7  [m_1]  •*  B£  = .78  * 10"7  [m-1] 

c.  B>  » 1 * 10-7  [m-1]  •*  B"  = 2.6  * lo'7  [m-1] 

c c 


In  these  three  cases  the  angle  of  arrival  « dependency  of  the  distance  D is  approximated  by 

a[degrees] 

a1.  D = 106  . a [kn]  a in  degrees 

b'.  D *>  420  . a tkm]  I 

D 

c'.  D *>  -135.  a [km] 
c 

The  path  delay  computations  At  = Ttota]_  - T0 

for  both  mclels  are  shown  in  Figures  5-7  . 


where: 

- T = total  path  delay  time  from 

total 

transmitting  to  receivirg  an- 
tenna 

- T = delay  time  of  a ray  propaga- 

ting at  the  earth's  surface  with 
a refractive  index  at  earth 

n = 1.0002  and  B'  = -1/R 
o a 

Figure 

The  figures  give  indication  of  larger 
At  values  in  the  spherical  model.  By  way  of 


► Ax  [ns] 

Figure  6.  Path  delay  as  function  of  a,  for 

B^  = -.78  * 10"7  tra’1]  (solid  line),  and 
B^  = .78  * 10  7 tm  *]  (dotted  line) 


B^  = -5.0  * 10  7 [m_1]  (solid  line),  and 
B"  = -3.4  * 10  7 [m  *]  (dotted  line) 

a 


— - At  [ns] 


Figure  7.  Path  delay  as  function  of  a,  for 

B^  = 1.0  * 10  7 [m  *]  (solid  line),  and 
B^  = 2.6  * 1C'7  [m-1]  (dotted  line) 


example,  differences  in  At  up  to  2.5  ns  are  found  for  angles  less  than  0.3  degrees  and  distances  less  than 
40  km. 
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To  verify  the  convergence  in  the  accuracy  of  the  numerical  computations  for  the  spherical  model  Ax 
computations  have  been  made  for  specific  angles  and  the  number  of  segments  (segment  thicknesses  which  have 
been  chosen  are  10,  5,  1,  and  0.25  meters)  has  been  increased.  In  Figures  8 and  9 the  accuracies  are  given 


elevation  angle  [degrees] 


elevation  angle  [degrees] 


e(x)  [ns]  - — 


— ► e (x)  [ns] 


Figure  8.  Error  e(T)  in  At  computations  as  function  Figure  9.  Error  e(T)  in  At  computations  as  function 


of  segment  thickness  for 
B'  =-5.0  * 10‘7  [.if1] 


of  segment  thickness  ds  for 
B*  = 1.0  * lO-7  [m-1] 


-7  _7  -i 

10  [m  ] and  B'  = 10  [m  ] . The  point  of  convergence  in  the  origin  has  been  obtained  by 


for  B’  =-5 

extrapolation  of  the  curve  for  a segment  thickness  of  1 m and  a thickness  of  0.25  m.  From  these  results  it 
can  be  seen  that  for  segment  thicknesses  less  than  10  meters  the  accuracy  lies  within  1 ns,  which  is  less 
than  the  differences  between  the  spherical  and  the  planar  model. 


4.  The  modified  logarithmic  refractive  index  profile 


An  application  of  the  theory  is  given  for  a refractive  index  profile  characteristic  of  a duct  above 
water.  This  model  has  been  chosen  because  strong  gradients  and  thus  strong  convergence  or  divergence  of  the 
rays  can  be  expected  just  above  the  water  surface.  At  the  same  time  there  is  the  possibility  to  inves- 
tigate whether  a simple  relationship  exists  between  the  amplitude  of  incoming  rays  and  their  delay  time. 

For  the  microwave  link  the  following  is  selected: 


3d  H 


2d  i 


- distance  between  transmitter  and  receiver  = 40  km 

- antenna  heights  = 80  m 


The  modified  logarithmic  refractive  index  profile 
above  a flattened  earth  is  given  by  (Figure  10). 


M(z)  = C, 


Jz  - Zj  - (d 

m * » 


+ zo>  • 


In 


where 


d = duct  height  (5-50  meters) 


- M(z) 


zq  = roughness  parameter  at  the  earth's  surface 
(10~3ra> 


Figure  10.  Modified  logarithmic  refractive  index 
profile  above  a planar  earth 


Zj  = height  where  M(Zj)  = rfj 


and  characterizes  dM/dz  at  high  altitudes 


(z  > 3d). 


20  30  40  50 

— ► <3[m] 


10  20  30  <0  50 

d[m] 

(a) 


10  20  30  40  50  0 10 

— - d[m] 

(b)  (c) 


iputational  results  as  function  of  duct  height  d. 

lei  2:  k = 4/3,  Mj  = 325;  a)  time  delay  difference  St;  b)  relative  amplitude  of  the  indirect 
r;  c)  angles  of  arrival 
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The  refractive  index  n above  the  spherical  earth  is  determined  by: 


M(z)  = (n  - 1 + —)  106  -»  dM/dz  = (dn/dz  + 1/R  ) 106 
R a 

a 


where  dn/dz  = refractive  index  gradient  above  a spherical  earth  with  earth  radius  R . 

a 


For  a fictive  earth  radius  R = kR  it  jo  found  that: 

f a 


f (dn/dz  + 1/R  ) ■*  dM/dz  = kR  ~ C2 
a a 


R_  = 


for  z > 3d 


-2  - 
6x  [ns] ' 


m — ' — i — 1 — i- > — i — 1 — r 


0 10  20  30  <0  50 

— - d[m] 
(a) 


(b) 


(c) 


Figure  13.  Computational  results  as  function  of  duct  height  d; 

Model  3:  k = 2,  Mj  = 325:  a)  time  delay  difference  6t:_  b)  relative  amplitude  of  the  indirect 
ray;  c)  angles  of  arrival 


This  means  that  k determines  the  modelconditions  for  high  altitudes.  The  computations  with  1 meter  segmen- 
tation and  for  k = 1,  4/3  (standard  atmosphere  above  the  duct)  and  k = 2,  are  presented  in  Figures  11,  12 
and  13.  In  these  figures  6t  is  the  delay  time  difference  and  r is  the  relative  amplitude  factor  of  the  two 
rays  in  this  two-way  propagation  model. 

A dominant  k-dependency  on  <5t  and  a clear  correlation  between  6t  and  ct  are  found  as  a function  of  the 
duct  height,  while  no  simple  relationship  exists  between  6t  and  r.  From  Section  3 it  is  known  that  the 
appropriate  spherical  model  has  to  be  used  for  accurate  5t  computations. 


Conclusions 


In  this  paper  it  is  proposed  that  for  path  delay  computations  in  multipath  fading  conditions  sphe- 
rical models  be  used.  In  the  spherical  model  analytical  formulas  for  the  ray  parameters  have  been  used, 
assuming  constant  refractive  index  gradients  per  segment.  By  doing  so  the  computing  time  can  be  reduced 
considerably  relative  to  numerical  solutions  of  the  differential  equations.  The  accuracy  analysis  shows 
strong  convergence  in  path  delay  results  by  decreasing  segment  thicknesses.  The  example  of  a two-way  sphe- 
rical model  based  on  surface  duct  layers  above  water  indicates  a strong  relationship  between  time  delay 
differences  and  angles  of  arrival  combined  with  a dominant  dependency  of  the  model  on  the  time  delay 
computations. 
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DISCUSSION 


M. P.M.  Hall  (U.K.):  Your  Figures  1-4  show  a rather  simple  case  where  both  transmitter  and  receiver  are 

situated  at  a height  above  which  the  refractive  index  decreases  sharply  with  height,  and  below  which  it 
increases  sharply  with  height.  Could  you  please  comment  on  how  the  main  conclusions  of  your  paper  would 
be  different  for  a more  general  ducting  geometry  in  which  the  terminal  height  would  not  be  at  this  region 
of  refractive  index  change,  and  where  the  curvature  of  each  ray  would  change  significantly  with  height  and 
even  change  direction  along  its  path? 

L.P.  Ligthart  (The  Netherlands):  In  the  spherical  model,  arbitrary  refractive  index  profiles  may  be  used 

and  antenna  heights  may  vary.  Maximum  time  differences  A t occur  in  the  case  you  mentioned  and  are  shown 
in  Figures  1-4.  The  main  conclusion  remains,  as  can  be  seen  from  Figures  11-13  in  which  a modified 
logarithmic  refractive  index  profile  has-been  taken. 

U.  Buse  (Germany):  At  which  frequency  range  are  your  models  valid? 

L.P. Ligthart  (The  Netherlands):  At  4 and  6 GHz. 

S.  Rotheram  (U.K.):  The  use  of  a flat  earth  with  an  earth  flattening  transformation  and  a modified  refrac- 

tive .index  should  give  the  same  answer  as  a spherical  earth  and  the  ordinary  refractive  index.  In  two 
dimensions  (cylindrical  geometry)  the  earth  flattening  transformation  can  be  made  exact. 

L.P.  Ligthart  (The  Netherlands):  The  earth  flattening  transformation  yields  a modified  refractive  index, 

which  gives  approximate  results.  In  the  limit  of  small  elevation  angles,  exact  results  are  obtained.  In 
comparison  between  planer  and  spherical  propagation  models,  differences  have  been  found  only  in  the  path 
delay. 

L.,  Boithias  (France):  Vous  ne  tenez  pas  compte  des  variations  horizontales  de  l'indice  de  refraction. 

Les  erreurs  introduces  ainsi  ne  sont-elles  pas  supSrieures  5 celles  que  vous  envisagez  de  corriger? 

L.P.  ligthart  (The  Netherlands):  Horizontal  variations  in  the  refractive  index  can  influence  the  numerical 

results.  Therefore  experimental  verification  of  the  models  is  needed. 

T. A.Th.  Spoelstra  (The  Netherlands):  This  is  a comment.  The  necessity  to  use  a planar  or  spherical  earth 

approximation  in  the  model  depends  on  the  accuracy  one  needs.  In  radio  astronomy,  correlation  of  signals 
from  telescopes  at  distances  of  a few  hundred  meters  make  a spherical  approximation  mandatory.  Furthermore 
it  should  be  valid  for  all  elevation  angles. 

N.  Amitay  (U.S.):  Roughly,  what  improvement  can  you  expect  to  obtain  by  using  angular  diversity? 

L.P.  Ligthart  (The  Netherlands):  Deep  fadings  can  be  reduced  between  5 to  10  dB  by  using  angular  diversity 

The  effects  on  group  delay  depends  on  delay  time  differences. 


